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ABSTRACT: We explore the sign problem in strongly coupled lattice QED with one flavor of Wilson 
fermions in four dimensions using the fermion bag formulation. We construct rules to compute 
the weight of a fermion bag and show that even though the fermions are confined into bosons, 
fermion bags with negative weights do exist. By classifying fermion bags as either simple or 
complex, we find numerical evidence that complex bags with positive and negative weights come 
with almost equal probabilities and this leads to a severe sign problem. On the other hand simple 
bags mostly have a positive weight. Since the complex bags almost cancel each other, we suggest 
that eliminating them from the partition function may be a good approximation. This modified 
partition function suffers only from a mild sign problem. We also find a simpler model which does 
not suffer from any sign problem and may still be a good approximation at small and intermediate 
values of the hopping parameter. We also prove that when the hopping parameter is strictly infinite 
all fermion bags are non-negative. 
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1. Introduction 

Strongly correlated many body fermion problems is an exciting area of research today [jl|, The 
main theoretical challenge in the field is to compute physical quantities of interest from first prin- 
ciples. Most methods that are currently used involve approximations that can be justified only in 
some regions of the parameter space. For problems where none of these approximations can be 
justified, the computational challenge is daunting. The Monte Carlo method is the only method 
which may be reliable in such cases. Unfortunately, this method also suffers from sign problems 
that arise due to the quantum nature of the underlying system ^ The final answer usually 
depends on delicate cancellations between many different quantum amplitudes which the Monte 
Carlo approach is unable to accomplish efficiently. The physics of nuclear matter and strongly 
correlated electronic systems are classic examples where the sign problem has hindered progress. 
Attempts to circumvent or solve the sign problem continues to be an important area of research and 
is also the focus of the current work. 

While a general solution to sign problems may not exist [j|], solutions have been found in 
specific cases when problems are reformulated using new variables. For example, while bosonic 
quantum field theories with a non-zero chemical potential suffer from a sign problem in the con- 
ventional formulation in the world line approach these sign problems disappear Even in 
fermionic quantum field theories, where the origin of the sign problem is the Pauli principle, new 
solutions are beginning to appear. In the conventional approach fermions are integrated out and the 
partition function is written in terms of bosonic degrees of freedom with a Boltzmann weight equal 
to the determinant of a matrix [g]. If this determinant is non-negative then the sign problem is ab- 
sent and today such problems can be solved using the popular hybrid Monte algorithm jl0| ] and its 
variants [pj]. On the other hand in many interesting cases the determinant can be negative or even 



complex. In such cases the conventional approach offers little hope for further progress. Recent 
research has shown that the world line formulations offer an alternative approach. Instead of inte- 
grating out the fermions at the beginning, considering their world lines and then re summing over 
only a limited class of these configurations leads to new solutions of the sign problems [|12|, |13| ]. 
The idea of using the world line aproach in two dimensional lattice field theories which usually do 



not suffer from sign problems has a long history [ 14, O, 16, 17, Nq]. Recently these developments 



have been unified under the framework called the "fermion bag" approach which shows that the 
new solutions to fermion sign problems can emerge in any dimension [19]. Basically one identifies 
independent dynamical regions over which the fermions naturally hop. These dynamical regions, 
called fermion bags, behave like non-local degrees of freedom. The weight of a fermion bag is 
just the path integral inside the fermion bag. When field theories are written in terms of fermion 
bags, sign problems may be absent since the weight of the fermion bags can be non-negative. The 
fermion bag approach allow us to solve some problems that seemed difficult or impossible in the 
conventional approach [|19[|, thanks to new algorithms |2C]. Being a relatively new idea not many 
examples have been studied and more work is necessary to understand the potential of the method. 

In this work we construct the fermion bag approach to four dimensional lattice QED with one 
flavor of Wilson fermions at strong gauge couplings. In a sense this is an extension of previous work 
in two [14] and three dimensions [18]. Wilson fermions contain a parameter called the hopping 
parameter referred to here as n. It is well known that the determinant of the one flavor Wilson Dirac 
operator in the background of a strongly fluctuating gauge field configuration can be negative for 
some values of k. Hence the conventional approach suffers from a sign problem in this region. 
Recently it was shown that the sign problem is absent in three dimensions when the partition 
function is written in terms of fermion bags [18]. Is this true in four dimensions? The current work 
was motivated by this question. At strong gauge couplings fermions are confined into bosons and 
fermion bags are regions where these bosons hop around. The weight of the bag is then a sum over 
all paths the fermions can explore within the bag while remaining confined. Since the fermions 
are always paired there is a possibility that the bags will have a non-negative weight. However, we 
show here that this is not the case. Fermion bags with negative weight do exist, suggesting that the 
underlying bosonic model remains frustrated. 

Although the fermion bag approach does not solve the sign problem, we can learn about the 
nature of the sign problem and some practical solutions from it. First, we can analytically prove 
that fermion bags with non-negative weights only contribute at k = oo. Thus, the fermion bag 
approach is able to solve the sign problem at this special point, while the conventional approach 
has a very severe sign problem there. Second, we find that at small n most bags that contribute have 
a positive weight. Negative weight bags begin to enter the partition function only for k > k c as 
in the conventional approach. Third, we find that large bags which are topologically simple (to be 
explained later) are also almost always positive. Large complex bags on the other hand have both 
positive and negative weights with almost equal probability. This creates a severe sign problem if 
they are allowed in the partition function. However, to a good approximation they seem to cancel 
each other and the partition function function is dominated only by simple bags. If one assumes this 
reasoning to be correct one obtains a new model that seems to capture at least some the interesting 
physics of the original model. This method of identifying new models by focusing on a class of 
fermion bags which capture important physics while being practically solvable may turn out to be 
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one of the main advantages of the fermion bag approach. 

Our paper is organized as follows. In section || we briefly review the sign problem in strongly 
coupled lattice QED with one flavor of Wilson fermions in the conventional approach. In section || 
we develop the fermion bag approach and construct diagrammatic rules to compute the weight of 
a fermion bag. In section ^ we classify bags as simple and complex and compute the weights of 
some small bags. We give examples of bags with negative weights. We also find the distribution 
of simple and some complex bags and use it to justify that complex bags do not contribute to 
the partition function. In section |5] we contruct a model without a sign problem that most likely 
contains the physics of parity breaking. We also give an analytic proof that the weight of fermion 
bags at k = oo are non-negative. Section ^ contains our conclusions. 



2. Sign Problem in the Determinant Approach 

Let us briefly review the sign problem in the conventional approach to strongly coupled lattice QED 
with one flavor of Wilson fermions. The partition function is given by 



[dipdip}[d(j)}exp(-S[^,ip,(t>}) (2.1) 
where the Wilson fermion action is given by 

S = ~Y1 + ^ x+a T a _e~^^ x ) + ij^Mfe (2.2) 

x,a x 

with the definition T± = (1 ± 7 a )/2. We denote the four Hermitian Dirac matrices as j a , a = 
1, 2, 3, 4. We also define 75 = —71727374 for later convenience. For explicit calculations we will 
use the chiral representation in which 

-(::)-(::)■ 

The four 2x2 matrices r a are defined by (ia, I) in the four vector notation. Note that a are the 
three Pauli matrices. The lattice fields ip x and t/j x represent the two independent Grassmann valued 
four component Dirac spinors on each hyper-cubic lattice site x and 4>x,a is the compact £7(1) 
lattice gauge field. In this work we choose open boundary conditions for convenience. Further note 
that our definition of k is two times the conventional definition of k [21]. 

The conventional approach is to integrate out the fermions and express the partition function 
as simply an integral over gauge fields. In this approach the Boltzmann weight of each gauge 
field configuration is simply the fermion determinant of the Wilson Dirac operator Dy/ [4>] in the 
background of that gauge field. More explicitly. 

Z = j [#] Det^LV [</>])• (2-4) 



where 



(D w [<t>]) x , y = Y1 S x+a>y T%e i ^ + S^y+aT^e-^ + ^6 x>y (2.5) 
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The Wilson Dirac operator satisfies the relation D^ w ^ = ^D w which can be used to show that 
eigenvalues of D w are either real or come in complex conjugate pairs. For n < 0.25 all real 
eigenvalues can be shown to be positive. However, for larger values of kappa there can in principle 
be an odd number of negative eigenvalues. Hence the determinant can be negative. The negative 
determinant is necessary to violate the Vafa-Witten theorem [^2|] and allow for the spontaneously 
breaking of the parity symmetry that is expected to occur for k > k c [23]. One expects k c ~ 0.5 
at strong couplings. 




Figure 1: Average value of the sign of Det(D w ) on 4 4 and 6 4 lattices as a function of k obtained using 1000 
random gauge field configurations. The mild sign problem on 4 4 lattice for 0.5 > k > 1.4 is just a finite 
size effect. 



In order to show the sign problem we compute the sign of the determinant of D\\/ on 4 4 and 
6 4 lattices in the background of a 1000 random U(l) gauge field. We plot the average value of this 
sign as a function of k in figure [T]. As expected the determinant approach encounters a severe sign 
problem when k > k c ~ 0.5. The sign problem continues to be severe even at k = oo. In this 
work we construct the fermion bag approach to this problem. 



3. Fermion Bag Approach 

At strong coupling we can first perform the link integral over the gauge field connecting x and 
x + a exactly to obtain an expansion of the partition function in terms of powers of Grassmann 
variables on each bond. We get 



/ 



%M^ne'U, + . + ? I+ „r 2 e-^) (3.!, 



We can also expand the exponential of the mass term on each site in terms of powers of Grassmann 
variables 

e -^ /K = yfiy[-^M\ (32) 

n=0 
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Collecting all the Grassmann variables on each site and performing the integration over the Grass- 
mann variables using the identity 



[dip] [dtp] (VOii^ji Ok rfj 3 0/0*3 ^hWu^. 



.74 



(3.3) 



we can rewrite the partition function as a sum over bond variables k x>a = 0, 1,2,3,4 and site 
variables n x = 0, 1, 2, 3, 4. We will refer to n x as the number of monomers on the site x and k x>a 
as the number of dimers on the bond connecting the site x and x + a. Thus, in the monomer, dimer 
representation the partition function given by 



Z= E H("B[n,k]) 

[n,k] B 



(3.4) 



where ojb is the weight of a "fermion bag" B which is simply the set of sites connected by k Xj0l ^ 0. 
Note that every site belongs to a unique bag and fermions only hop within the sites of the bag. The 
Boltzmann weight of a fermion bag, ojb, is the sum over a well defined set of fermion hoppings 
within the bag. Of course there is no need for lob to be positive. However, in this work we prove 
that lob is indeed positive when k = and |«| = oo. We also find evidence that to a good 
approximation certain class of fermion bags, which almost always have positive weights, dominate 
the partition function for a range of values of k. 

Let us now construct the rules for calculating lob- For this purpose we define 
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it is easy to show that Y°L = S\ a S + , a , r° = S_ a S-, a for every a. Further 



I ( (I - sis 2 T ai rl 2 ) 



R SlS2 




(3.6) 



(I— siS2T ai T a2 )/2 is a 2 x 2 matrix which can be parametrized as F^ s ^ 2 exp(in* 1 i ; 
(t(7t/4)). table [I] lists the values of F and n for various possibilities. 



where R*^ a2 



2 . 

a 2 





«2 


01,02 


CI1,CI2 


a 
4 

i = 1,2,3 
i = 1,2,3 


a 

i = 1,2,3 
4 

i=l,2,3 


5(1 - S!S 2 ) 

1 

V2 
1 

V2 
1 

V2 




nu = sis 2 5 ik 
nk = -sis 2 5ik 
n>k = —s\s 2 €ijk 



Table 1: Values of and n* 1 *^ that enter the definition of R s a \% ■ 



Note that R^ a = and R a a = I, while for all other values of a\ and a 2 , the matrix Rf£ 



S1S2 

0.2 



is (1/ y/2) times a (1/2, 0) representation of an 0(4) rotation matrix. Using these relations we can 
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write 



k,l 

[Yl ( S +^3 m (^+«)m (ipx+Jn (Sl )Ct )r 



I 

m,n 



(3.7) 



The integration over the Grassmann variable then leads to specific rules that help to compute the 
weight Uflofa bag. 

The weight of a bag turns out to be the trace of the product of Dirac tensors associated to each 
site. The explicit form of these Dirac tensors are discussed below. But first is useful to remember 
the constraint that every site in the bag must satisfy 

n x + ^2 k x>a + k x ,-a = 4 - (3-8) 
a 

Here we have defined k x - a = k x - a ^ a . Based on the allowed values of n x , each site in the bag can 
be one of seven types as shown in table |2[ We call these as type-0,l,2,3,4a,4b and 4c depending on 
the number of dimers attached to the site. Note that there are three types of sites with four dimers 
attached to it. We distinguish them because the rules to compute the weights are slightly different 
for each of them. We also use two types of diagrammatic representation for each vertex: a detailed 
diagram and a minimal diagram. The detailed diagram shows each fermion line and is helpful in 
the actual computation, while the minimal diagram just shows the dimers (or a single monomer 
when no dimers exist on the site). Given the minimal diagram the detailed diagram can be uniquely 
obtained. 

The simplest site is type-0 site where n x = 4. Such a site forms its own bag since it is not 
connected to any dimers. It has a weight 

u B = W = k" 4 . (3.9a) 

Next consider the type-1 site with n x = 3 and k x ^ siai = 1 where si = ±1 and a\ is one of four 
possible positive directions. The contribution to the weight of the fermion bag due to such a site is 
in the form of a Dirac tensor (Wi)?*"* and is given by 

K -3 



in 3 



S- sl , ai St ua \ . =0. (3.9b) 



Hence a fermion bag cannot contain a type-1 vertex. 

Next consider the type-2 vertex where n x = 2 and k x>siai = k XjS2a2 = 1. In this case the 
Dirac tensor associated with this site is of the form (^^yjjj" 2 that contributes to ujb is given 
by 

0^^)i 1> i 2 ;j 1 j 2 — n\~(^—si,ai )iifci {S-S2,a2)i2k2 £ kik2ksk4 £ hl2k:jk4 (S sit ai )hji (si 

= K (S- Sljai S s a ) ( 5'_ S2jQ , 2 S'^ 1 ) = K (-R Q , 1 1 a 2 )jij 2 (-^a2ai)*2il 
V / tiJ2 \ ' «2J1 

(3.9c) 
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Note that if a\ = «2 and s\ = S2 the tensor is zero. 

Next consider the type-3 site with n x = 1 and k XjSiai = k X:S2a2 = k XjS , sas = 1. Then 




— IK ^(^ Q ,\a 2 )*li2(-^0203)*2i3(^«3«l)*3il 

+ (^Q 1 id 3 )h 33 (Ra 3 a 2 )*3J2 )*2 jl J (3.9d) 

Note that again all the dimers must be in different directions otherwise the site weight is zero. 

Finally we can have a type-4 site with n x = 0. In this case we have four directions given by 
k x ,siai = k X)S2Ct2 = k X)Sm3 = k XjSia4 = 1. Now there are three possibilities: Type-4a is one in 
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which all the four dimers are in different directions and then we get 



m/ a, \ siai ' S2a2 ' S3Q;3 ' S4Q4 — (<3 \ (Q \ fc \ fq \ e 

\ VV 4 /ii,i2»*3,»4lJl)j2j3>j4 ~~ ^ 8 U a lHlkl \ J -S2,Ot 2 )l2k2 \°-S 3 ,Ot 3 )l 3 k3\ lJ -S4,a4jl4k4 t -klk2k3k4 

£ hhhh i^li ,ai)hh (^S2,a 2 )hj2 {^s 3 ,a3 )hj3 {^34^4)1434 
(B SlS2 )■ ■ (E S2Sl )■ ■ (E S3S4 V • (E S4S3 V • 

_(DS lS2 \. . (r>S 2 S 3 \. . (r,S 3 S4 \. . (JDS4S1 \. . 
1 (RS1S3 y . ( ns 3 si y . /ns 2 s 4 y . I raS4S2 y . 

_(jysiS3 y . fos 3 s4 y . (n>s 4 s 2 y . fpsasi y . 

V JL ai«3^ij3 V-"'a3a4^3j4 V-"'a4a2'' J 4j2 V- 1 L a 2 oi / J 2jl 
_ ( 'osis 3 y . (0S3S2 y . (f>S2S4 y . ( r>S4Si y . 

_(T>S1S4 y . (OS4.S3 y . / Z? S 3 S 2 ). . (ftS- 



3S2S1 
l a2ai 

(3.9e) 

Type-4b is the site where k X;Siai = 2 and k XjS3Ce3 = k X;S4a4 = 1. The above expression then 
simplifies to 



f'T/f/ b "\ SlQl ' S3a3 ' S4Q4 — _ 1 (R SlS 3 V • ('p s i s 4 y . _ (j? s i s 3 y . (T>siS4 y . 

\ 4 -'*l,i2,»3)»4;ilj2iJ3j4 2 V Q l Q 3^1-?3 l >- lt Oi04'' 4 2j4 V-' l «i«3 A2J3 V- ,t Oi0 4 Mj?4^ 

C7?*3Si y . /oS4si y . _ C D«3Si y . / 0S4S1 y . ^ 

V""'a3ai/«3J1 V- ,t 040iy«4J2 I -"'0301 7*3.72 V- ,t 04ai 7*4.71 J 

(3.9f) 

Here the extra factor of 1/2 is due to the fact that there are two dimers on one of the bonds and 



this leads to the an extra factor l/(2!) 2 present in Eq.(3T ). This extra factor can be divided equally 
between the two vertices that the dimer connects. Type4-c site is obtained if k XjSlCei = k XjSsa3 = 2. 
In this case we get 

(^)S J1J2J3J4 = \{K\ S a\)\r 2 ) lll2 (r 2 ) jli2 (r 2 ) i3i4 (r 2 ) i 3 i4 (3.9g) 

Again the extra factor of 1/4 is due to two double dimers and the factor Fa\fa 3 comes from the R 
terms. This completes the classification of all the vertices 



4. Sign Problem with Fermion Bags 

Using the rules of the previous section it is possible to compute the weights of fermion bags nu- 
merically. However, it is exponentially difficult to compute the weight when the bag contains many 
type-3 and type-4a sites. In order to make progress, we label the bags with the number of type-3 
and type-4a vertices it contains. Thus, a bag of type (7x3, 77.4) contains 773 type-3 sites and 714 type- 
4a sites. The (0, 0) bags contain no type-3 and type-4a sites and will be referred to as simple bags. 
Bags in which either 77,3 or 77,4 is non-zero will be called complex bags. Below we will argue that 
this classification in terms of (773, 774) helps in understanding the origin of the sign problem. 
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Bag 




dimer representation 


bag type 


[ 




□ 


] 


1 —4 
64 K 


(0,0, 0,0; 1,1) (0,0, 0,0; 2,1) (1, 0, 0, 0; 1, 2) 
(1,0,0,0;2,1) (0,1, 0,0; 1,1) (2, 0, 0, 0; 2, 2) 
(1,1,0,0;1,2) 


(0,0) 




















(0,0, 0,0; 1,1) (0,0, 0,0; 2,1) (1, 0, 0, 0; 1, 1) 
(1,0, 0,0; 2,1) (0,1, 0,0; 1,1) (0, 1, 0, 0; 2, 1) 












21 -12 


(4,1) 










128^ 


(2, 0,0,0; 2,1) (1,1,0,0;1,1) (1, 1, 0, 0; 2, 1) 
(0,2, 0,0; 1,1) (2,1,0,0;2,1) (1, 2, 0, 0; 1, 1) 
























(0,0, 0,0; 1,1) (0,0, 0,0; 2,1) (0, 0, 0, 0; 3, 1) 
(1,0, 0,0; 2,1) (1,0, 0,0; 3,1) (0, 1, 0, 0; 1, 1) 
(0,1, 0,0; 3,1) (0,0, 1,0; 1,1) (0, 0, 1, 0; 2, 1) 
(1,1,0,0; 3,1) (1,0,1,0;2,1) (0, 1, 1, 0; 1, 1) 












3 -8 
32^ 


(8,0) 
















_5_ K -14 
2048 ^ 


(0,0,0,0;l,l) (0,0, 0,0; 2,1) (0, 0, 0, 0; 4, 1) 
(1,0, 0,0; 2,1) (0,1, 0,0; 1,1) (0, 1, 0, 0; 3, 1) 
(0,1, 0,0; 4,1) (0,0, 0,1; 2,1) (1, 0, 3, 0; 2, 1) 
(1,0, 3,0; 3,1) (1,1,0, 0;3,1) (0, 1, 1, 0; 1, 1) 
(1,1,3,0; 3,1) 


(2,2) 



Table 3: Some small fermion bags and their weights. 



By now it should be clear that every fermion bag can be uniquely represented through the 
dimers of the bag. We represent these dimers using the notation (x\, x%, X3, X4; a, k) where Xj 
represent the four dimensional coordinates of the site inside the bag from which k dimers emerge 
in the positive direction a. Some examples of fermion bags, their dimer representation and their 
weights are given in table ^| Although all the bags shown in the table have a positive weight we do 
find bags that have both zero weight and negative weights. However these bags are more complex. 
Two examples of negative weight bags and one example of zero weight bag are given in table |] 
along with their weights : Bag-1 is a simple bag which contains twelve type-2 and two type-4b 
vertices. Bag-2 is a complex bag made up one type-4a, seventeen type-2, four type-3 vertices. 
Bag-3 is a simple loop bag with zero weight. 

In order to understand the sign problem we have generated fermion bags of a fixed type at 
random on an L 4 lattice using a worm algorithm. In the case of simple bags, we exclude single site 
bags and plaquette bags for convenience. We then analyze the probability distribution of bags of a 
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Ba, 


?1 






(0,0,0,0;1, 


1) 


(0,0, 0,0; 2,1) 


(1,0, 0,0; 2,1) 


(1,0, 0,0; 3, 2) (0,1, 0,0; 1,1) 


(1,0, 1,0; 4,1) 


(0,0, 1,0; 1, 


1) 


(0,0, 1,0; 4,1) 


(1,0, 1,1; 4, 1) 


(0,0, 1,1; 2,1) 


(0,1, 1,1; 3,1) 


(0,0, 1,2; 1,1) 


(0,0, 1,2; 3, 


1) 


(0,0,2, 1; 2,1) 


(0,0,2, 1; 4,1) 














Ba, 


?2 






(0, 0, 0, 0; 2, 


1) 


(0,0, 0,0; 3,1) 


(0,1, 0,0; 2,1) 


(0,1, 0,0; 3,1) 


(0,0, 1,0; 1,1) 


(0,0, 1,0; 2,1) 


(0,2, 0,0; 3, 


1) 


(0,1, 1,0; 1,1) 


(0, 1, 1,0; 2, 1) 


(1,1, 1,0; 2,1) 


(1,0, 0,0; 2,1) 


(1,0, 0,0; 3, 1) 


/ 1 i n n. a 

(1, 1,0,0; 4, 


1) 


(1,2,0,0; 1, 1) 


z' i on n. o 1 A 

(1,2, 0,0; 3, 1) 


(1,1, 0,1; 4,1) 


(2, 2, 0,0; 3,1) 


(1, 1,0,2; 1, 1) 


(1,0,0,2; 1, 


1) 


(1,0, 0,2; 2,1) 


(2, 1,1,0; 2,1) 


(2, 1,1,0; 4,1) 


(2, 0,0, 2; 2,1) 


(2, 1,0,1; 3,1) 


(2, 1,0,1; 4, 


1) 
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(1,0, 1,0; 4, 


1) 


(0,0, 1,0; 1,1) 


(0,0, 1,0; 4,1) 


(1,0, 1,1; 4,1) 


(0,0, 1,1; 2,1) 


(0,1, 1,1; 3,1) 


(0,0,1,2;1, 


1) 


(0,0, 1,2; 3,1) 


(0,0, 2,1; 2, 1) 


(0,0, 2,1; 4,1) 







Table 4: Examples of bags with negative weight (Bag 1, Bag 2) and zero weight (Bag 3). Bag 1 has a 
weight of -1.220703125000 x 10~ 4 and Bag 2 has a weight -1.430511474609 x 10~ 6 at n = 1. 



given type using the bag action density defined by 

*B = — -^log(l"B|), (4-1) 

where Ng is the number of sites in the bag. In figure ^ we plot the distribution of (0, 0), (2, 0) and 
(2, 1) bags on the 2 4 lattice with open boundary conditions as a function of Sb- For each type of 
bag we have generated 10 4 bags. The left panel contains the distribution of positive weight bags 
while the right panel shows the distribution of the negative weight bags. We find that all simple 
bags (or (0, 0) type bags) turn out to have positive weights. On the other hand complex bags ((2, 0) 
and (2, 1) type bags) do contain negative weight bags. In the (2, 0) case we find 6987 positive 
and 3013 negative weight bags, while in the (2, 1) case we find 5983 positive and 4017 negative 
weight bags. We have repeated a similar analysis on a 5 4 lattice where we have generated more that 
3 x 10 4 bags. These results are plotted in figure [|. In this case a small number of simple bags do 
have negative weights. But the positive and negative weight complex bag distributions are almost 
identical for both (2, 0) and (2, 1) bags as can be seen from the figure. 

Based on figures ||] and ^ we conclude that as n-^ and increase (in other words as the bags 
become more complex) the distribution of positive and negative weight bags become more and 
more identical and hence the sign problem becomes severe. On the other hand simple bags are 
dominated by positive weight bags. Thus, we believe that to a very good approximation complex 
bags will cancel each other and the partition function is dominated by simple bags. Assuming this 
to be true an interesting effective model of strongly coupled QED emerges in which the partition 
function only contains simple bags. This model may share some of the physics of the original 
model. On the other hand it may be studied in its own right since it will have a much milder sign 
problem. We postpone this study to a future publication. 
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Figure 2: Distribution positive weight bags (left panel) and negative weight bags (right panel) as a function 
of the action density Sb on a 2 4 hyper-cubic lattice with open boundary conditions. Three types of bags are 
shown: (0, 0)-type (top) (2, 0)-type (center) and (2, l)-type (bottom). See text for more details. 
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Figure 3: Distribution positive weight bags (left panel) and negative weight bags (right panel) as a function 
of the action density Sb on a 5 4 hyper-cubic lattice with open boundary conditions. Three types of bags are 
shown: (0, 0)-type (top) (2, 0)-type (center) and (2, l)-type (bottom). See text for more details. 



5. Fermion Bags with Non-negative Weights 

Can we construct a model of strongly coupled QED with Wilson fermions which is completely 
free of the sign problem in the fermion bag approach? In order to answer this question we identify 
fermion bags with non-negative weights. We find that there are three classes of fermion bags for 
which we can prove analytically that the Boltzmann weights are always non-negative. The first is 
the trivial bag consisting of a single site for which the weight is Wq = k~ 4 . The second class 
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are loop bags which have a loop topology. These bags only contain sites of type-2. Since they are 
closed loops of confined fermion and anti-fermion world lines, their weight is a square of a trace of 
an SU(2) matrix and hence real and non-negative. Interestingly if we modify the original action to 

5 = "E fe^in^ 1 '"^^ + ^ x+a i> x+a ^ x+a r«e- i ^ x ) + -^ x ^ x (5.i) 



x.a 



it is easy to argue that only trivial bags and loop bags are produced and the sign problem is com- 
pletely solved. Since at small values of k type-3 and type-4 sites are naturally suppressed, this 
model may be a good approximation to the original model at small and intermediate values of k. 



It most-likely contains the parity breaking phase transition of the original model Q23Q. The hand- 
waving argument is as follows: At small values of k the loops are small while at large values loops 
proliferate the entire lattice and hence are naturally large. It is easy to show that on a finite lattice 
(ip^ystp) = due to the parity symmetry. On the other hand the two point correlation function 
(^xl^x ifryl^y) will be non-zero. This two point correlation function gets contribution through 
an open loop with the end points at x and y. Intuitively, at small values of k, since the loops 
will be small, the correlation function decays exponentially to zero for large separations. On the 
other hand at large values of k, when the loops are large, the correlation function will decay as 
a power law and thus signaling the spontaneous breaking of parity. This phase transition can be 
studied efficiently using a worm-type algorithm. We postpone this study to the future. It would be 
interesting to understand the nature of this transition. Note that the above model will suffer from 
a severe sign problem in the conventional approach since additional auxiliary fields in addition to 
the usual gauge field will have to be introduced to convert the action into a fermion bi-linear. This 
is yet another example of a model which is solvable in the fermion-bag approach rather than the 
conventional approach. 

The third class of bags with non-negative weights consist only of type-4 sites. These bags 
arise naturally when k = oo. The proof that the Boltzmann weight is non-negative is a bit more 
involved and relies on the bi-partite nature of the lattice. Let us briefly sketch the proof here. From 



Eq.(3.9e) we know that the contribution to the weight from each site within the bag comes from 



the tensor 

h,i2,i3M'Jlj2,j3,ji (^— Sl,Ql)llfcl ("S'-S2,02)«2fc2 W~ S3,(X3 ) 13^3 (^—34 ,04 )l4k4^ki k 2 k 3 k4 



(Ssi,ai)j 1 l 1 {Ss2,a2)j2h (Ss3,a3)j 3 l 3 {$84,014) jihEhhhh 



(5.2) 



If we define 



rpS\ai,S2CH2,S30t3,S40t4 ( q \ (q \ I q \ Iq \ (S.'W 

11,12,13,14 ~ \ -si,ai)iiki\ —S2,a2)i2k2\°-s 3 ,a 3 )i 3 k3\iJ—S4,a4)i4k4 t kik2k3k4 V J -~V 

we see that 

yi/" _ rpsiai, 3202,3303,3404 (rpSiai,S2a 2 ,S3a3,S4a4\* (5 4") 

4 il,*2,«3,M \ jl,J2J3,J4 ' \ J - ) 

This structure of W4 shows that, on a bi-partite lattice, the Boltzmann weight of the bag will be the 
square of the magnitude of a complex number obtained by tracing over the product of T's on each 
site. 

Although the above argument proves that all the fermion bags with type-4 vertices will have 
non-negative weights, as far as we know, a practical Monte Carlo algorithm seems impossible due 
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to the fact that it will be exponentially difficult to compute the Boltzmann weight of large fermion 
bags. In a sense, the sign problem may still be hidden in this computational difficulty. 



6. Conclusions 

In this work we have constructed the fermion bag approach to strongly coupled lattice QED with 
one flavor of Wilson fermions in four dimensions. We found that at k = oo all fermion bags have 
non-negative weights. On the other hand fermion bags with negative weights do exist and create 
a severe sign problem at intermediate values of k. By classifying bags as simple and complex 
we could show that complex bags almost cancel each other in the partition function while simple 
bags are almost always positive and hence contribute to the partition function. This suggests a 
simple solution to the sign problem. We simply approximate the partition function as the sum of 
contributions form simple bags. 

This approximate solution to the sign problem is similar in spirit to the meron cluster approach 



[13]. There special clusters called meron clusters appeared with equal weight but opposite sign in 
the partition function. Allowing meron clusters in the partition function would create a very severe 
sign problem. However, since they come with exactly equal weight and opposite signs, they cancel 
exactly and thus the sign problem was solved completely. In the current situation, the cancellation 
of complex bags is only approximate and suggestive. So, while we cannot justify rigorously that it 
is correct to ignore them in the partition function we believe it to be correct. In the future it would 
be interesting to study the partition function generated by simple bags alone. 



Finally, we have also constructed a simpler model (Eq. (p. 1|)) that consists of loop bag and 
does not suffer from the sign problem. Simple arguments suggest that this model contains two 
phases : a parity symmetric phase at small values of k and a phase where parity is spontaneously 
broken at large values of k. It would be interesting to study the nature of this phase transition in 
three dimensions. 
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A. Fermion Determinant versus Fermion Bags 

In this appendix we provide checks that confirm the correctness of the rules that we constructed 
in section ||, to compute the weights of fermion bags. We compute the partition function on a 
small lattice by integrating Det(D w ) over the gauge fields exactly and identifying contributions 
from each of the fermion bags that are produced in the process. For simplicity we choose a 3 x 2 
and a 3 x 3 lattice on the xy — plane as shown in figure |J The fermion bags that are produced 
in these small lattices already capture all types of vertices of table ||. We label each lattice site 
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by an index % and assign different values k = K{ to each site. The partition function will then 
be polynomial in Kj's such that powers of K\ are related to the number of monomers on site i. 
This helps identify terms in the partition function as weights of the fermion bags. Finally we set 
Kq = K\ = K2 = ■ ■ ■ = k to compute the weight of the bag. 

8 7 6 



A 








1 





Figure 4: Lattices on which we compute the partition function exactly. The site labels are used to label 
different values of k which help to identify the contribution to a particular fermion bag. 

On a 3 x 2 lattice we get find the partition function to be a sum of 1 1 terms given by 
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Collecting the terms into 8 categories each of which contributes to the weight of a fermion bag, we 
can compute these bag weights. In the top eight rows of table || we compare the weights computed 
by this method with the one computed using the fermion bag rules of section ||. In this calculation 
we find all vertices listed in table ^| except for type-4a vertex. To find a bag with type-4a vertex 
we compute the partition function on a 3 x 3 lattice. In order to simplify the calculation, we set 
the links between sites 2 and 1, 2 and 3, 8 and 5, 8 and 7 to be zero. Then the partition function 
contains 6 terms and are given by 



2564 

^3 ^6 ^7 ^8 



+ 



Ki-^ Ki^K/^ A^g fcrj /^g 



+ 



4t--2 „ 2 ^.4 ^.4 ,-2 ^.4 ^.4 ^.4 2 
012 3 4 6 7 8 5 
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/l i, ■ 4 4 j.. 4 i, . 2 j.. 2 i, . 2 i, . 2 ,,. 4 4 O f-i 4 1,-4 1,-4 1,-4 1,^ 4 

4: Aq A A^2 A ^ A,j rign^^n/g 5 ZjOUAq A^ 2 8 5 



+ 



444444444 

01 2 3 4 ^ Ay Ag A 



(A.2) 



These can be divided into 4 categories, 3 of them contribute to bags that have already been enu- 
merated within the top eight rows of table ||| The only bag which includes type-4a vertex is given 
in the last row of table Bl These results confirm the rules constructed in section B. 
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Table 5: Comparison between the determinant approach and the bag approach on a 3 x 2 lattice (top 8 rows) 
and on a 3 x 3 lattice (bottom row) with open boundary conditions. Each diagram corresponds a unique term 
in the partition function. 
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